OPEN 9 ACCESS 



C 



International 

Journal 

Of Modern Engineering Research (IJMER) 

) 



Mhd Flow of A Non-Newtonian Fluid Through A Circular Tube 

NASER R. NASER \ R. K .KHARE 2 ’ AJIT PAUL 3 

*(. Department Of Mathematics and Statistics Shiats Allahabad India ) 

**( Department Of Mathematics and Statistics Shiats Allahabad India ) 

***(Department Of Mathematics And Statistics Shiats Allahabad India ) 



( \ 

ABSTRACT : The work deals with the study of unsteady motion of a dusty fluid through a 

uniform pipe with sector of a circle as cross-section, the fluid and the dust particle clouds are 
static at t = 0. the influence of time of the flow is dependent to the pressure gradient. We get the 
solution by analysis of the partial differential equation by variable separable method and Laplace 
transform. Finally with the aid of graph, the results had been analyzed 
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I. INTRODUCTION 

The MHD problem becomes more important rapidly during the past few years, it is used in fluidization, 
use of dust in gas cooling systems for heat transfer process, in environmental pollution, in dust collection, in 
nuclear reactor cooling, etc. P.G Saffman[12] has discussed the stability of the laminar flow of a dusty gas in 
which the dust particles are uniformly distributed. Liu[13] has studied the flow induced by an oscillating infinite 
flat plate in a dusty gas. N.C. Ghosh [3] obtained the solution of MHD flow of a visco-elastic fluid through 
porous medium. P. Samba siva Rao [13] has got the solution for unsteady flow a dusty viscous liquid through 
circular cylinder. This paper deals with the study of laminar flow of an unsteady dusty fluid through pipe with 
sector of circle as cross-section.P.D.Verma and A.K.Mathur[14] studied the unsteady flow of a dusty viscous 
liquid in a circular tube.Bikash Chandra Ghosh and N.C.Ghosh[15] obtained the solution for MHD flow of a 
visco-elastic fluid through porous medium.M.H.Hamdan[16] has obtained the dusty gas flow through porous 
medium. The equation of motion of an unsteady viscous incompressible fluid with uniform distribution of dust 
particles are given by P. G. Saffman (1962). 



II. FORMULATION AND SOLUTION OF THE PROBLEM 

V.u= 0 (continuity) (1) 
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V.t?=0 (continuity) (3) 

— + (- 5 . V’)t? =— (u — v) (Linear Momentum) (4) 

dt ■ - m 



Where U and V are the velocity of fluid and dust , t is the time,p is density of the gas , p is pressure of the 
fluid, N is the number density of dust particles, V is kinematic viscosity, k is Stoke’ s resistance is mass of 
the dust particle, is the magnetic field parameter , is the conductivity, C is non-Newtonian factor . 

We assume an unsteady laminar flow of an incompressible viscous fluid with uniform distribution of dust 
particles through porous medium in a uniform pipe with sector of a circle of radius CL as cross-section .The flow 
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is due to the influence of time dependent pressure gradient. Both the fluid and the dust particle clouds are 
supposed to be static at the beginning .It is assumed that the dust particles are spherical in shape and uniform 
in size and the number density of the dust particles is taken as a constant throughout the flow .The velocity 
distributions of the fluid and dust particles are given by 

«r = 0 , u e = 0 , u s = u z (r, 0, t) 

v r =0, v s = 0 , v z = v z (r, 6, t) (5) 




Figure. 1: Geometry of the flow in circular channel 



By equations (5) , (2) and (4) we get 
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Consider the following non-dimensional quantities 

u s — au, v s — av , p — pvp, z — az , 

„ _ - . _ k — mv 

r — ar . t — — r . k — —r 
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The initial conditions are 



it(r,0,z) = 0 ,v(r, 9, z) = 0 when t = 0 ,r G (0,1) and 9 E (—a, a ) 
(9) 

The no slip boundary conditions are 
t fc (0,O3) 

u{r, 9, z) = 0 ,v(r, 9,z) = 0 When r = 0 and r = 1 
«(r, 6, z) = 0 ,t>(r, 0, z) = 0 When 0 — ±a 



Let 



( 6 ) 
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( 8 ) 



( 10 ) 
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We give the definition for Laplace transformations as follows 

u — e~ xt u dt and v — e~ xt v dt 

Applying the Laplace transformation to equations (7) and (8) and 
using boundary conditions, we get 



xu=p(X)+(— + — +—) + -(v-u) - 



vp Cl+c a ) 



U 



(ID 



( 12 ) 



XV = \ (u — v) 
0 



(13) 



The boundary condition now becomes 



u = v = 0 When r = (0,1} and 0 = [—a, a} 
Where p(A T ) is Laplace transform of 3 (O' 
Equation (13) implies 

_ U 

v = 

l-¥x 

By (12) and (15) Then 

d‘u . 1 3 u . 1 d‘u -, 2 - . _ , x _ 
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e assume the solution of equation (16) in the following form 
u{r, 8 ) = w^r, 8 ) + w 2 (r) 

Substitution of u(r, 0) in equation (16) 

S 2 v 3“v>\ 1 / 3w ± dw.A 1 3 V, 2 1 i % 

-^T + 77f + + -^ ) + + W Z ) 

+P(x)= o 
If w 2 satisfies then 

“ M 2 w 7 + P(x) = 0 

or r or 



And wq satisfies 

or m r g r -r“ 0 &* 

If u(r, 8) is inserted in no slip boundary conditions, someone can obtain 

u(0, 0) = w 1 (0, 0) + w 2 (0) = 0 , tt(l, 8) = w^l, 8 ) T w 2 (l) = 0 , 
u(r,a) = w^T'Ct) +w 2 (r) = 0 , u(r f —a) = w-^r,— a) -fw 2 (r) = 0. 

The equation (19) is the Bessel function whose solution is 
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w 2 (r) = C^iMr) + C 2 K Q (Mr) + 



P(x) 



Af* 



( 21 ) 



Where J 0 , K a are the zeroth order modified Bessel function of first and second kind 
respectively and C v C 2 are constants. 

By conditions w 2 (0) = 0 and w 2 (1) = 0 in equation (21) 

Where ( u r ,u 9 ,u s ) and (v r „v e ,v s ) are velocity components of fluid and dust particles 
respectively. Above will be shown in the figure. 1 



We get 



P(x) 


A (Afr) 
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( 22 ) 



By this equation we see that the velocity at the centre of the cylinder is zero. Using variable 
separable method, the solution of the problem (20) With conditions 
WjXO, 9) = 0 ,w ± (l ,9) - 0 r w^(r r a) = -w 2(r) 

^(r, —a)= —w 2 (r) is obtained in the form 

= (AI x (Mr) 4- BK x (Mr))(C cos(A0) + D sin(^e)) (23) 

Where l x and K x are the modified Bessel-function of order 2. Of first and second kind 
respectively and A,B,C and D are constants by the same way B = 0 since the velocity at the 
centre of the cylinder is zero , K x (Mr) — > co as r — > o 
Then the solution is given by 
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Substituting (22) 
£t(r,0) — 



and (24) in (17) we obtain 
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Using u in equation (15) can see that 
v(r, 9 ) - 



P(x) 


r^Ofr)-^Cw)l 
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(24) 



(25) 



(26) 



Case of Impulsive motion 

Let 3(0 — Where S(t) is Dirac delta function and p 0 is constant. 



Now equation (25) and 
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By taking the inverse Laplace transform of u and v then 



u 



= y fj 
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e * ir (U-a: L ff) 0) 

i+tl+a^) 3 iHKl+* a £) a 



Where (m= 1,2,3, ) are the positive zeros of / 0 («) = 0 

Now r = 0,1 and 9 = +a respectively are 
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(30) 



Where 

Xl 2/>(l+ C a )v£ ’* 2 2/j(l+C a )v£ 

flj - P(1 + C 2 ) V (1 + l - J?Q + 

R 2 - i 2 - 4p(l + C 2 )vJ?(a^ z a 2 -p(l + C z )va m 2 ) 

Table No.( 1) 

(T=l , /?=(). 05 , m=2 , Po=0.6 , p=2 , c=3 , cr=l , a=l) 
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s. 

No. 




Fluid 

velocity 

t=0.2 


Fluid 

velocity 

t=0.4 


Fluid velocity 
t-0.6 


Fluid velocity 
t=0.8 


1 


0 


0.5314 


0.5879 


0.6505 


0.7198 


2 


0.5 


0.5298 


0.5845 


0.6451 


0.7121 


3 


1 


0.5247 


0.5746 


0.6293 


0.6893 


4 


1.5 


0.5165 


0.5584 


0.6038 


0.6530 


5 


2 


0.5052 


0.5365 


0.5700 


0.6055 


6 


2.5 


0.4911 


0.5098 


0.5294 


0.5497 


7 


3 


0.4745 


0.4791 


0.4839 


0.4888 


8 


3.1 


0.4709 


0.4725 


0.4744 


0.4762 


9 


3.2 


0.4672 


0.4659 


0.4648 


0.4637 




We see in above table (1.1) the fluid velocity increase as the magnetic field increase , for every parameter (t) 
.the reason of this increasing is the value of terms X lf X 2 which are power for the constant (e) in the 
numerator then the value for the fluid velocity will be increase 



I IJMER I ISSN: 2249-6645 I 



www.ijmer.com 



I Vol. 5 I Iss. 8 I August 2015 I 64 I 



Preparation of Papers for International Journal of Modern Engineering Research 



Table No. (2) 

(T=l , /? = 0.01 , m=2 , Po=0.6 ,p= 2 , c=3 , cr=l , a=l) 



S. No. 




Dust velocity 
(=0.2 


Dust 

velocity 

at=0.4 


Dust 

velocity 

t=0.6 


Dust 

velocity 

t=0.8 


1 


0 


0.5180 


0.5734 


0.6345 


0.7021 


2 


0.5 


0.5167 


0.5704 


0.6296 


0.6949 


3 


1 


0.5127 


0.5618 


0.6153 


0.6739 


4 


1.5 


0.5063 


0.5477 


0.5922 


0.6404 


5 


2 


0.4974 


0.5285 


0.5615 


0.5965 


6 


2.5 


0.4862 


0.5050 


0.5244 


0.5446 


7 


3 


0.4730 


0.4779 


0.4827 


0.4875 


8 


3.1 


0.4701 


0.4721 


0.4739 


0.4758 


9 


3.2 


0.4671 


0.4662 


0.4650 


0.4639 



OO 

CZ> 

CO 

CO 

cz> 

c\T 

cz> 

II 




magnetic field (BO) 



Similarly in this phase we see in above table(1.2) the dust velocity increaseas the magnetic field increase, for 
every parameter (t) .the reason of this increasing is the value of terms X lf X 2 which are power for the 
constant ( e) in the numerator then the value for the dust velocity will be increase 
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Table No. (3) 

( p a = 0.6 , 0=0.05 , m=2 ,1=1 ,<7=1 , p= 2 ,a=l) 



S.No. 


Non-newtonian 


Fluid velocity 


Fluid velocity 


Fluid velocity 


Fluid velocity 




factor (c) 


t=0.2 


t=0.4 


t=0.6 


t=0.8 


1 


0 


0.0891 


0.0261 


0.0076 


0.0022 


2 


0.1 


0.0904 


0.0267 


0.0079 


0.0023 


3 


0.2 


0.0945 


0.0288 


0.0088 


0.0027 


4 


0.3 


0.1013 


0.0324 


0.0104 


0.0033 


5 


0.4 


0.1107 


0.0376 


0.0128 


0.0043 


6 


0.5 


0.1224 


0.0446 


0.0162 


0.0059 


7 


0.6 


0.1363 


0.0535 


0.0210 


0.0083 


8 


0.7 


0.1518 


0.0644 


0.0273 


0.0116 


9 


0.8 


0.1687 


0.0772 


0.0354 


0.0162 


10 


0.9 


0.1865 


0.0918 


0.0452 


0.0223 


11 


1 


0.2047 


0.1080 


0.0570 


0.0301 




Table (3) reveals as the non-newtonian factor(c) of fluid increases, the velocity of fluid phase 
increases for every time parameter (t) , because (c) is in the denominator of x lr x 2 in (29). 
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Table No.( 4) 



S.No. 


Non-newtonian 
factor (c) 


Dust velocity 
t=0.2 


Dust velocity 
t=0.4 


Dust velocity 
t=0.6 


Dust velocity 
t=0.8 


1 


0 


0.1224 


0.0376 


0.0110 


0.0032 


2 


0.1 


0.1299 


0.0385 


0.0114 


0.0034 


3 


0.2 


0.1344 


0.0410 


0.0125 


0.0038 


4 


0.3 


0.1417 


0.0453 


0.0145 


0.0046 


5 


0.4 


0.1515 


0.0515 


0.0175 


0.0059 


6 


0.5 


0.1637 


0.0596 


0.0217 


0.0079 


7 


0.6 


0.1777 


0.0698 


0.0274 


0.0108 


8 


0.7 


0.1931 


0.0820 


0.0348 


0.0148 


9 


0.8 


0.2095 


0.0960 


0.0440 


0.0201 


10 


0.9 


0.2264 


0.1116 


0.0550 


0.0271 


11 


1 


0.2434 


0.1285 


0.0687 


0.0358 



(P 0=0.6 , /?=0.05 , m=2 , 1=1 , cx=l ,p = 2 , a=l) 




Table(4) reveals as the non-Newtonian factor(c) of dust increases, the velocity of dust phase increases for 
every time parameter (t) , because (c) is in the denominator of ^1^2 in (30).The graph (1) , (2) ,(3) and(1.4) 
supported the previous tables 
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III. RESULT AND DISCUSSION 

From the first and second table, we can show that the relation between the fluid velocity and 
magnetic field is Inverse relation. Similarly for the dust, from the third and fourth table we can show the relation 
between the fluid velocity and non-newtonian factor is proportional relation .similarly for the dust. By previous 
we can show the relation between magnetic field and non-newtonian factor is inverse. 
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